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Single slepton production in association with a top quark at the CERN Large Hadron 
Collider (LHC) is one of the important processes in probing the R-parity violation cou- 
plings. We calculate the QCD next-to-leading order (NLO) corrections to the pp —> 
tl~(t£ + ) + X process at the LHC and discuss the impacts of the QCD corrections on 
kinematic distributions. We investigate the dependence of the leading order (LO) and the 
NLO QCD corrected integrated cross section on the factorization/renormalization energy 
scale, slepton, stop-quark and gluino masses. We find that the uncertainty of the LO cross 
section due to the energy scale is obviously improved by the NLO QCD corrections, and 
the exclusive jet event selection scheme keeps the convergence of the perturbative series 
better than the inclusive scheme. The results show that the polarization asymmetry of 
the top-quark will be reduced by the NLO QCD corrections, and the QCD corrections 
generally increase with the increment of the t\ or g mass value. 
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I. Introduction 



Supersymmetry (SUSY) is one of the most appealing theories as an extension of the stan- 
dard model (SM). It provides an elegant way to solve the gauge hierarchy problem and cancel 
the quadratic divergences of the radiative correction to the Higgs boson mass. The mini- 
mal supersymmetric standard model (MSSM) contains SM particles, their superpartners and 
an additional Higgs doublet. In order to avoid the rapid proton decay, a discrete .R-parity 
symmetry 

HQ 

is introduced in the MSSM, which implies a conserved quantum number, 
R p = (_i ) 3B + L + 25, ; with B, L and S being baryon number, lepton number and spin of the 
particle, respectively. For all the SM particles R p is equal to +1, while for all the superpartners 
R p is — 1. So this symmetry leads to the result that the superpartners can only be produced in 
a pair and the lightest SUSY particle (LSP) is stable. However, such a stringent symmetry ap- 
pears to be a theoretical basis; especially, we know that a stable proton can survive by imposing 
either L- or ^-conservation 3]. At the same time, there is not enough experimental evidence 
for .R-parity conservation, so i?-parity conservation is not necessary in the MSSM. Moreover, 
non-zero -R-parity violating (RPV) couplings could provide small neutrino masses, which could 
explain the phenomena of neutrino oscillation experiments. Thus, there are strong theoretical 
and phenomenological motivations to introduce partial i?-parity violation into the most general 
representation of superpotential with the respect to the renormalizability and preservation of 
the gauge symmetries of SM and supersymmetry, which can be written as 4|-[2] 

l - A A A A 1 „ A A ,s A 

= -e ab X ijk L^L b j E k + eabXijkLIQ^Dk + -e a p 1 \ ijk U^D j D\ + e ab ^iL1Hl, (1.1) 

where i, j, k denote generation indices, a, b (= 1,2) are SU(2) isospin indices, and a, /?, 7 
are SU(3) color indices. Li (Qi) are the left-handed leptons (quarks) SU(2)-doublet chiral 
superfields, and Ei (Ui, L>i) are the right-handed leptons (up- and down-type quarks) SU(2)- 
singlet chiral superfields. H\ i2 are the Higgs chiral superfields. The Xijk, A^ jfc and /ij are 
the dimensionless L-violating coupling coefficients, and X'-j k are the i?-violating dimensionless 
coupling constants. As mentioned above, B and L cannot be violated at the same time. In this 
paper we concentrate ourselves only on the L-violating couplings, where the coefficients \'^ k 



may be assumed to be non-zero. 

The LHC can be used as a top-factory, and it is advantageous to study the production 
and decay of top quark. We suspect that the processes related with top quark can be used to 
probe the new physics effects since the top mass is close to the weak scale [7J, particularly in 
the single top-quark production process the chiral structure of the interaction of the top quark 
may impact the polarization observable of the final top quark 8[. Therefore, the predictions 
including higher order corrections to single top-quark production within and beyond the SM 
are very important in exploring the new physics. 

The single charged slepton production in association with a top quark at hadron collider 
is induced by the non-zero term X' i3k LiQ 3 D k in Eq. dl.ip . The pp — > t£^(i£f) + X process at 
the LHC receives the contributions from the partonic processes gdk(gdk) — > t^(t^t)j where i 
and k are the generation indices. The second term in Eq. (II. lj) is related to the Born process 
pp — > t£^(i£^) + X, which can be expressed as 



C 



LQD 



ijk 



ViidkRdjL + djidkRViL + d* kR Vi L dji 



C-ibdkRUjL — UjidkR^-iL — d* kR l C iL UjL 



+ h.c. 



;i-2) 



For the NLO QCD calculations, the interaction vertices of two squarks and a slepton may be 



involved, which can be extracted from the general soft SUSY-breaking Lagrangian 



9] 



rsoft 



Kjk A \ y i d-iL dt R - hi u iL dt R ) + h.c 



i a jL a kR l Li UjL a kR 



;i.3) 



In the above equation it is assumed that the soft breaking terms have a universal dimensionful 
parameter A and are proportional to the dimensionless coupling constant A^ fc . In this work we 
take the SUSY-breaking parameter A = 0, ±1 TeV, separately. 

The single top-quark pro duct ionpro cesses at colliders in the i?-parity violating MSSM has 



been studied in several references 
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1CH16I|. The top-quark production in association with 



HE? 



a slepton £j at hadron collider has been studied at leading order (LO) in Refs. [17j,ll8(, there 
the authors performed the signal analysis, and found that the final states in the t£~(i£f) 
production at the LHC have distinct kinematic signatures, which can be distinguished from the 
backgrounds. 
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In this paper, we present the calculations of the next-to-leading order (NLO) QCD correc- 
tions to this process. The paper is organized as follows: In section II, we present the calculations 
for the relevant partonic processes and parent process pp —> t£j~(i£f) + X at the LO and QCD 
NLO. In section III, we give some numerical results and discussions. Finally, a short summary 
is given. 

II. Calculations 

A. LO calculation 

In both the LO and NLO calculations, we apply FeynArts3.4 and FormCalc5.3 packages 
Q Q to generate F_ d lagra ,„s, their corresponding a m p.itudes, and to — the 
amplitudes, separately. In Table [1] the upper 2a bounds on X' i3k originating from Refs. [18|, |21[ 
are listed. There the coefficients A^ 33 < O(10~ A ) stem from assuming m u < leV and left-right 
mixing in the sbottom sector. Since we have the strong constraints on the A^ 33 coupling shown 
in Table [1] and the low (anti)bottom luminosity in parton distribution function (PDF) of the 
proton, which indicates there cannot be any significant production rates via gb(gb) — > t£~(t£'l) 
partonic processes at the LHC, we ignore their contributions in the following calculation. 



^131 


0.019 x (m ir JlOOGeV) 


y 

A 132 


0.28 x (m ir JlOOGeV) 


y 

A 231 


0.18 x (m ir /100GeV) 


y 

A 232 


0.45 (m~ SR = WOGeV) 


A 331 


0.45 (rriq = lOOGeV) 


y 

A 332 


0.45 {nig = lOOGeV) 


^.33 


O(10~ 4 ) 







Table 1: Upper 2cr bounds on A^ 3fe , where m,- I(B) is the mass of the left (right) handed squark 
Ql(r) ■ 



Due to the CP-conservation the production cross section for gdk — > t£~ (i — 1, 2, 3, k — 1, 2) 
subprocess is the same as that for the gdk — > tt[ [i = 1, 2, 3, k = 1, 2) subprocess. In this 
section we present only the calculations of the former subprocess. There are two tree-level 
Feynman diagrams contributing to the partonic process of gdk — > tt^ (i — 1, 2, 3, k — 1, 2) as 
shown in FigfjTa) (for the s-channel) and FigJIJb) (for the t-channel). 

The expression of the LO cross section for the partonic process g{pi)dk{p2) — > t{pz)^ ip^) 
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(a) (b) 
Figure 1: The tree- level Feynman diagrams of the gd k — >• t£~ partonic process. 



has the form 



i i (2-jr) 4 : r color 

a (gd k ^t£r) = ^—2- I J2\M LO (gd k ^ t£r)\"dn 2 , (2.1) 

spin 

where the factors \ and ^ come from the averaging over the spins and colors of the initial 
partons, respectively, s is the partonic center-of-mass energy squared, M.Lo{gdk t£^) is the 
amplitude of the tree- level Feynman diagrams shown in FigJTJ The summation in Eq.f l2.ip is 
taken over the spins and colors of all the relevant initial and final particles. The phase-space 
element dQ 2 is expressed as 



(2.2) 



The LO cross-section for the parent process pp — > gd k — > + X at the LHC can be 
obtained by performing the following integrations: 



a LO = ^2 / dxi / dx 2 a (gd k -> U i ) [G g/Pl (xi, fif)G dk /p 2 (x 2 , fif) + (x x <H> x 2 ,P\ P 2 )] , 

d k =d,s J ° Jo 

(2.3) 



where Gj/a(x, /i/) is the PDF of parton j (= g or d k ) in proton A (= Pi,P 2 ) which describes 
the probability to find a parton j with momentum xPa in proton A, pLf is the factorization 
energy scale. We adopt the CTEQ6L1 PDFs in the LO calculations. 
B. Real and virtual corrections 

In the NLO calculations we use the dimensional regularization method in D = 4 — 2e 
dimensions to isolate the UV and IR singularities. Some of the virtual QCD one-loop diagrams 
are shown in Figf2J 




Figure 2: Some of the QCD one-loop Feynman diagrams for the gdk — > t£^ partonic process. 
The upper indices x,w — 1,2, and the lower indexes i and k run from the first generation to 
the third generation. 

The NLO QCD corrections to the gd k —¥ t£~ partonic process can be divided into two 
components: the virtual correction and the real radiation correction. There exist ultraviolet 
(UV) and infrared (IR) singularities in the virtual correction and only IR singularity in the real 
radiation correction. The IR singularity includes soft divergence and collinear divergence. In 
the virtual correction component, the UV divergence vanishes by performing renormalization 
procedure, and the soft IR divergency can be completely eliminated by adding the contribution 
from the real gluon emission partonic process gdk — > t£i + g. The collinear divergence in the 
virtual correction can be partially canceled by the collinear divergences in real gluon/light- 
(anti) quark emission processes, but there still exists residual collinear divergence which will be 
absorbed by the redefinitions of the PDFs. 

The one-loop diagrams can be divided into two independent parts. One is the SM-like 
component arising from the diagrams including gluon/quark loop, another is the pure super- 
symmetric (pSUSY) QCD part where each diagram includes a gluino/squark loop. Correspond- 
ingly we divide the counterterms also into the SM-like QCD and the pSUSY QCD parts. The 
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definitions of the relevant counterterms are adopted as 



Ql/r 



1 + l-SZ 



1 + 2 { SZ L/R 



q,(SM-like) 



+ 5Z 



q,(pSUSY)\ 
L/R J 



Ql/r 



. r . x (SM-like) . x (pSUSY) 

m t — y m t + om t —m t + om t + om t 



_ ^ Z (SM-like) _|_ SZ (pSUSY)^ 



g s ^ 9s + Sg s = g s + 



(SM-like) 



+ Sg. 



(pSUSY) 



(2.4) 



where g s is the strong coupling constant, qi/R denote the fields of top and <4(= cf, s) quarks, 
and represents the gluon field. We renormalize the relevant fields and top-quark mass in the 
on-shell scheme {22]. For the renormalization of the QCD strong coupling constant g s and the 



i?-parity violating coupling coefficient A- 3fc , we use the MS scheme 
of the g s can be expressed as 

(SM-like) 

-ln- 

g s Aix 2 e 3 



23 



24( | . The counterterm 



Sg 



9^ 

(pSUSY) 



Air 
a s (n r 



g (SM-Uke) 1 1 2 1 2 

2 e 3 3 fxf. 



i=1.2 



i=i,2 



— + — ln-| + > — In— ^ + > — In — 

2 e 3 /i r 2 12 „V 12 ^ 

U=u,c,t D=d,s,b 



2 1 



(2.5) 



where 1/e = l/e^-7B + M47r), /?( 5M -^ e ) = | nj/ - | an d /3( pSC/sy ) = _§j\r_ I( n , / + 2 ) 

with AT = 3 and riif = A. 

With the Lagrangian shown in Eq.f ll.2p . the counterterm of the ti — d^—t vertex is expressed 
as below: 

(2.6) 



= -x*. ( ^ + **** + ^' 



a 2 



By using the MS" scheme to renormalize the £i — dk — t coupling, we get 

0A i3k — A i3k 1 2 / 7re' 



(2.7) 



where Cp = 4/3. After the renormalization we get a UV-finite virtual correction to the partonic 
process gdk —> ttT . 



7 



Figure 3: The tree-level Feynman diagrams for the real gluon emission process gdj, — > t£ i + g. 



The real radiation correction includes the contributions from the gluon and light- (anti)- 
quark emission processes. The contribution of real radiation processes is at the same a s order 
as the virtual correction to the partonic process adj. — > ti~ in perturbation theory according to 



the Kinoshita-Lee-Nauenberg(KLN) theorem 



25 



26| . The tree-level Feynman diagrams of the 
real gluon emission partonic process g{pi)dk{p2) — > t(p^)£^~(p/i)q(p^) are depicted in Figj3) We 



27] to isolate the IR singularities for 



adopt the two cutoff phase-space slicing (TCPSS) method 
the real emission subprocesses by introducing two cutoff parameters 5 S and 5 C . The arbitrary 
small soft cutoff 5 S separates the three-body final state phase space of real emission subprocess 
into two regions: the soft region (E§ < 5 s \/I/2) and the hard region (E s > 5 s V§/2). The 
collinear cutoff 5 C separates hard region into the hard collinear (HC) region and hard non- 



collinear (HC) region. The region for real hard gluon/light- ( ant i) quark emission with s 15 (or 
S25) < 8 c s (where = (pi + Pj) 2 ) is called the HC region. Otherwise it is called the HC 
region. Then the cross section of the real gluon emission partonic process can be written as 

af(gd k U i9 ) = a s g + of = of + of c + of (2.8) 

where Og, of c and of c are the cross sections in the soft gluon region, hard collinear region 
and hard non-collinear region, respectively. 

The light- (anti) quark emission contribution at the QCD NLO to the pp — > t£~ + X process 
includes the partonic channels: (1) q(q)dk — > + q(q), (2) gg — >■ tl~ + dk- The corresponding 
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Figure 4: The tree-level Feynman diagrams for the real light-(anti)quark emission processes. 
(1-8), (9-10) and (11-12) are for the gg — > tiid^, qd k — > t£~q and qd k — > t£^q (q = u,d,s,c) 
partonic processes, respectively. 

Feynman diagrams of these partonic processes at the tree-level are shown in FigJU There the 
diagrams by exchanging the identical incoming quarks in FigsJH(9)-(10) for the partonic process 
qdk — >• t£~q, i.e., q = d k {= d, s) are not drawn. 

Again we use the TCPSS method to splitting the three-body phase space into collinear (C) 
and non-collinear (C) regions, and get the cross sections for the subprocesses gg — > t£~d k and 
q(q)dk —> t£^q(q) at the tree level expressed as 

* R (gg^UA)=ti=ai + drg, (2.9) 
a R (q(q)d k ^ U^q)) = fifa = a% + a%. (2.10) 

The cross-sections in the non-collinear region, a£, a^-^in Eqs.f l2.9p . (I2.10p ) and a^ c (in 
Eq.f l2.8l) ). are finite and can be evaluated in four dimensions by using the general Monte Carlo 
method. After summing all the contributions mentioned above there still exists the remaineding 
collinear divergence, which will be absorbed by the redefinition of the PDFs at the NLO. 
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C. Total NLO QCD correction 

Given the NLO correction components to the cross sections of the subprocesses above, the 
full NLO QCD correction to the cross section for the pp — > ti^ + X process at the LHC is 
formally given by the QCD factorization formula as 

°nlo(pP -*ti~ + X) 

^ / dxi / dx 2 { [G j / Pl (x 1 ,fif)G k /p 2 (x 2 ,iJ,f)d-( jk )(x 1 X2s)] + (xi -H- x 2 ,Pi -H- P 2 )} , 
,,,, Jo Jo 



(2.11) 



where j and k sum over all possible types of initial partons contributing to the subprocesses 
up to the QCD NLO, i.e., (jk) represents (gd), (gs), (gg), (yd), (qs), (qd) and (qs) (where 



q = u,d,c,s). We consistently adopt the CTEQ6m PDFs|28|, |29] for Gj/p^p 2 )(xi, [if) and 
Gk/p 2 (Pi)( x 2, A*/)- The total NLO QCD correction Aa NLO to the pp — > tlj + X process can be 
divided into two-body term and three-body term, i.e., Aa^po — + Aa^ 3 \ The two-body 

term consists of the virtual correction and the cross section for the real gluon/light-(anti) quark 
emission process in the soft and hard collinear phase-space region. The three-body term consists 
of the cross section for the real gluon/light-(anti) quark emission process in the hard non- 
collinear phase space. Finally, the QCD corrected total cross section for the pp — > tl~[ + X 
process is 

onlo = ctlo + a (2) + a® . (2.12) 



III. Numerical results and discussion 

In this section we present and discuss the numerical results for the LO and NLO QCD corrected 
cross sections for the pp — > t£^(t£f) + X process at the early (y/s = 7 TeV) and future 
(y/s = 14 TeV) LHC. In order to check the correctness of our LO calculations, we list our 
LO numerical results for pp — > gd k — > tl~[ + X process in Table [2j and compare them with 



those presented in Table [2] and Fig. 3 of Ref. 



181 ] . Our results are obtained by employing the 



same input parameters and PDFs as used in previous work [18]. We can see that they are in 
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m h {GeV) 


pp -> UJ + X via A^ 31 = 0.1 


-> t£j + X via A = 0.1 


our a LO (fb) 


a LO (fb) in Ref. [18 


our cr LO (fb) 


olo (fb) in Ref.[18] 


200 


377.8(5) 


378 


73.5(4) 


73.5 


500 


51.80(7) 


52 


7.64(4) 


7.7 


800 


11.58(8) 


11.6 


1.370(7) 


1.4 



Table 2: The comparison of our LO numerical results for the pp — > gdk — > tij + X process at 
the y/s = 14 TeV LHC with those in Ref. 18j. The relevant parameters and the PDFs being 
the same as used in Ref. 18 . 



agreement. But we should say that most of the data taken from Ref.jlSJ] are read from the 
figure with large errors. 

In the following numerical calculation, we take one-loop running a s and two-loop running 



a s in the LO and NLO calculations, respectively 30|]. The number of active flavors is taken 
as N f = 4, and the QCD parameter Af° = 215 MeV and the CTEQ6L1 PDFs are adopted 



in the LO calculation, while A 



MS 



NLO calculation 



28 



326 MeV and the CTEQ6M PDFs are used for the 



29]. We set factorization scale and renormalization scale equal and take 
fi = jj,Q = (m t + m|.)/2 in default. The CKM matrix is set to the unit matrix. We ignore the 
masses of electron and u-, d-, s-, c-quarks, and take m t = 172.0 GeV and m h = 4.2 GeV in the 
numerical calculation 30]. 



In the calculation of the process pp — > ti~ + X we assume that the masses of the sleptons of 
three generations are degenerated with the values of = m^, (i = 1,2,3), and take the ^ p 
coupling parameters as (1) A- 31 = 0.1 and the other A' = 0; (2) A- 32 = 0.1 and the other A' = 0. 



We follow the SPA benchmark point SPSla' 



3l|, where the input parameters are m = 70 GeV, 



mi /2 = 250 GeV, A = -300 GeV, sign(n) = +1, tan/3 = 10, and m t = 172.0 GeV, but we 
assume that there is no left- and right-squark mixing in the first two generations and the 
degenerat sleptons having the masses mi = 189.9 GeV (i = 1, 2, 3). Then we get the left-right 



stop (sbottom) mixing angle 9 t = 57.026° (6b = 19.75°) by adopting the ISAJET program 
and the other relevant SUSY parameters with the values: 



32|, 



m-g = 607.1 GeV, m h = 366.5 GeV, m h = 585.5 GeV, 



m 



506.3 GeV, m~ b2 = 545.7 GeV, m Ul = m Sl = rriu R = 547.2 GeV, 
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ras, = ma 



564.7 GeV, mi = m Sl = 



546.9 GeV, 



m S2 = m^ L = 570.1 GeV, m^o 



97.7 GeV, fi = 396 GeV 



(3.1) 



The SM parameters used in the calculation are taken as follows |3Qf ] 



a, 



ew 



1/137.036, sin 2 9 W = 0.23119, m w = 80.398 GeV, 



m z = 91.1876 Gel/, m 



0.1057 GeV 



(3.2) 



In the following calculations we take the parameters stated above if there is no other statement. 

The independence of the full NLO QCD correction on the two cutoffs, 5 S and 5 C , is confirmed 
numerically. FigsJ5]^a,b) demonstrate that the total NLO QCD correction to the pp — > gd — > 
t£~ + X process with A = TeV at the y/s = 14 TeV LHC does not depend on the arbitrarily 
chosen values of the 5 S and 5 C within the calculation errors, where we take A^ 31 = 0.1 and other 
A' = 0. In Fig. 5 (a), the two-body correction Aa^ 2 \ three-body correction Aa^ and the total 
QCD correction (Aa^Lo) for the pp — > gd — )■ t£^+X process are depicted as functions of the soft 
cutoff 5 S running from 1 x 10~ 5 to 1 x 10 -3 with 5 C = 5 s /50. In Figj5]^b), the amplified curve for 
Acjnlo is depicted. The results in these two figures demonstrate that the total QCD correction 
to the pp gd — >• ti^ + X process is independent of 5 S and 5 C . It verifies the cancelation of 
the soft /collinear IR divergence in the total QCD correction to pp — >• gd — >■ tl~ + X. In further 
numerical calculations, we fix 5 S = 5 x 10 -5 and 5 C = S s /50. 

For the gluon/light-(anti) quark jet event selection we adopt two selection schemes: (1) 
Inclusive jet event selection scheme. With this scheme we accept all the real gluon/light- 
(anti)quark emission events. (2) Exclusive jet event selection scheme, by which we accept the 
real gluon/light-(anti) quark emission event satisfying the restriction of either p 3 ^ 1 < p^j et on 
the jet transverse momentum or rf et > r(?f t on the jet rapidity. In further calculations with 
exclusive jet event selection scheme, we take the cut parameters for gluon/light-(anti) quark jet 
as p C T*j et — 50 GeV and = 3. 

One of the main reasons to base the LHC analyses on higher order predictions is the stabi- 
lization of the dependence on the unphysical renormalization and the factorization scales. In the 
upper figures of FigsJ6](a,b,c), we show the dependence of the LO and NLO QCD corrected cross 



12 



(a) 



-Acj 

■Aa" 

-Ao 



-3000 - 
■6000 - 



Vs = 14TeV 
M~=189.9GeV 



-4.0 

Log(8J 



(b) 



2 710- 



\'s = 14TeV 
M- = 189.9GeV 



-4 

Log(5 s ) 



Figure 5: (a) The dependence of NLO QCD corrections to the pp — >■ gd — > + X process on 
the soft cutoff 5 S at the = 14 Tel/ LHC with A = 7W, = 189.9 GeV and 5 C = S s /50. 
(b) The NLO QCD corrected total cross section ctatlo with calculation errors as the function 
of S 8 . 

sections for the processes pp — >■ gd — )■ £i~ + X, — >• — >• t5£+ + X, and pp gs tt^ + X 
on the factorization/renormalization scale (ji/fio) with A = TeV at the >/s = 7 TeV^ and 
a/s = 14 TeV^ LHC, respectively. The corresponding K-factor [K = (Tnlo/glo] are shown in 
the lower figures of Figsj6]^a,b,c). Since the luminosities of the s- and s-quark in proton are the 
same, the observables for the process pp —> gs — )■ i£f + X are equal to those for the process 
pp — > gs — > tlj + X, we present only the plots for the process pp — >■ gs —> t£~ + X. There we 
assume \i = /i r — Hf for simplicity, and set A- 31 = 0.1 and other A' = in Figsj6|a,b), while 
\> i32 = 0.1 and other A' = in Fig^c). The curves labeled NLO(I) (NLO(II)) denotes the 
NLO QCD corrected cross section with the inclusive (exclusive) jet events selection scheme. 
The curves in the lower figures of Figsj6]^a,b,c) labeled (i), (ii), (iii) and (iv) correspond to 
the K-factors for (i) the inclusive scheme with yfs = 14 TeV, (ii) the inclusive scheme with 
y/s = 7 TeV, (iii) the exclusive scheme with a/s = 14 TeV, (iv) the exclusive scheme with 
yfs = 7 TeV, respectively. These notations are also adopted in the following figures. We can 
see from Figsj6]^a,b,c) when the scale \i runs from 0.1/xq to 4/i , the curves for the NLO QCD 
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corrected cross section becomes more stable in comparison with the corresponding curves for 
the LO. It demonstrates that the NLO QCD corrections can improve the scale uncertainty 
apparently, and the exclusive scheme keeps the convergence of the perturbative series better 
than the inclusive scheme in the plotted /x/z-to range. 

In FigsJ7]^a,b,c) we plot the LO, NLO QCD corrected cross sections and the corresponding K- 
factors as the functions of the final slepton mass rug. with A = TeV for the pp — >■ gd — >■ t£~+X, 
pp — >■ gd — >■ i£f + X, and pp —> gs t£~ + X processes, respectively. In each figure the curves 
for the a/s = 7 TeV and y/s = 14 TeV LHC are depicted. In Figsf7^a,b) we set the _R-parity 
violating coupling coefficients as A^ 31 = 0.1 and the other A' = 0, while in FigJT^c) we have 
X' i32 = 0.1 and the other A' = 0. We can see that both the LO and NLO QCD corrected cross 
sections decrease with the increment of the value of m~ t . . The curves in these figures show that 
the cross sections for t£~ and iif production processes are different. Unlike the luminosities 
of s- and s-quark in proton, being equal, the initial d-quark has a higher luminosity than the 
d-quark. This induces the cross-section of the process pp — >• gd — > t£^ + X to always be larger 
than the process pp — >■ gd — » t£f + X. 

The SUSY partners t\ and g appear only at one- loop level (shown in FigfS]). We show 
the virtual influences of and trig on the NLO QCD corrections by taking three different 
values of the SUSY-breaking parameter A (i.e., A = TeV and ±1 TeV) in FigsJSHTUl and 
FigsnU-EU separately. FigsErtMa,b,c) (FigsrjTHI^a,b,c)) show the NLO QCD corrections 
(Ao-jVLo) an d the corresponding K-factors as functions of {m g ) for the pp — >■ gd — )■ t£~ + X, 
pp -> gd -> t£f + X , and pp gs ^ t£j + X processes with A = TeV and ±1 TeV at the 
early and future LHC, respectively. There we adopt both the inclusive and exclusive schemes 
in the NLO calculations. In FigsJH]4T3]^a,b) we set A^ 31 = 0.1 and the other A' = 0, while in 
FigJ8HT3l(c) we take A^ 32 = 0.1 and the other A' = 0. In those figures most of the curves for the 
NLO QCD corrections and the corresponding K-factors always increase with the increment of 
m h ( m g) except those AaNLO curves in FigJTQTc) and the curves in the region of m g < 200 GeV 
in Figs fT^T a.b.c). It shows that the NLO SUSY QCD corrections to these processes generally 
increase in the large or m g region. That is because the SUSY QCD correction coming 
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Figure 6: The LO, NLO QCD corrected cross sections and the corresponding K-factors versus 
the factorization/renormalization scale (fi/no) with A = TeV at the LHC. (a) For the pp — > 
gd — > t£~ + X process, (b) For the pp — >■ gd — >■ + X process, (c) For the pp — )■ gs — >■ + X 
process. The K-factor curves labeled (i), (ii), (iii) and (iv) are for (i) the inclusive scheme with 
y/s = 14 TeV, (ii) the inclusive scheme with y/s = 7 TeV, (iii) the exclusive scheme with 
y/s = 14 TeV, (iv) the exclusive scheme with y/s = 7 TeV, respectively. 
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Figure 7: The LO, NLO QCD corrected cross sections and the corresponding K-factors versus 
the mass of slepton with A = TeV at the LHC (a) For the pp — >■ gd — )■ + X process, 
(b) For the pp gd t£f + X process, (c) For the —> gs — >■ + X process. The 
descriptions for the K-factor curves labeled (i), (ii), (iii) and (iv) are the same as in FigsJHJ 
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from the renormalization counterterm of the amplitude 8Ai involves the logarithm terms of 

(m? \ / m 2 \ 

-^r I and In f -j^-J at the SUSY QCD one-loop level. The logarithm term contributions 

can be absorbed by a redefinition of the coupling A' 33j. In this work we do not need to 

adopt this decoupling scheme, since in our consideration we take the squark and gluino mass 

quantitatively at the same order as the scale u,. 

Figs|S]and FigsJTTlshow that if we take A = 0, the NLO QCD corrections are still related to 
mi and m g . This can be understood from the fact that the NLO QCD contributions from the 
counterterms shown in Eqs. (j2.4j) — ( I2.5P and the loop diagrams such as Figsf2]^5,6), are relevant 
to the masses of t and g. We can also see the Aa NLO curves versus in Fig JTOT c) for the 
pp — > gs — » tlj + X process, and the Autvlo curves in the region around m g ~ 200 GeV 
in Figs JT5T a.b.c) are obviously distorted by the contributions from the diagrams involving the 
tj-dk—li coupling with A = 1 TeV. In order to understand the contributions from the non- 
zero ij—dif- li coupling more clearly, we provide the plots of 5a^ LO (= <tnlo{A = ±1 TeV) — 
&nlo(A = TeV)) versus the m g in Figs JT4T a.b.c) for the processes pp —¥ gd — » tlj + X, 
PP ~ >* dd ~ >* t^t + X and pp — > gs — >■ t£~ + X, separately. There we can see that in the region 
of m g < 400 GeV, 5a~^ LO is always positive, while Sa^ LO remains negative. 

In the following calculations we fix the SUSY-breaking parameter A = TeV. We show 
the polarization asymmetries (A t ) of the (anti)top quark as a function of m|. at the LO and 
NLO QCD for the processes pp — > gdk(gdk) t^Ti^t) + i n FigsJToTa.b.c.d). A t is defined 
as A t = , where N + and N_ refer to the numbers of positive and negative helicity 

(anti)top quarks respectively. In Ref.jl^], the LO A t at the LHC has been plotted and it shows 
that the LO A t changes signs for a slepton mass of around 870-900 GeV . As we can see in 
Figs JToT a.b.c.d). the polarization degree of the (anti)top quark has been reduced obviously by 
the NLO QCD correction. That is because the NLO QCD radiation corrections destroy the 
chiral structure of the interaction of top quark and reduce the LO polarization asymmetry A t . 
Although the cross sections for the processes pp — >■ gs — > i£f + X and pp — > gs — > tij + X are 
equal, those two mutually conjugate processes have opposite A t values as shown in Figs JToT c) 
and (d). 
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Figure 8: The NLO QCD correction to the LO cross sections (Acttvlo) and the corresponding 
K-factors versus with A = TeV at the \fs = 7 TeV and \fs = 14 TeV LHC (a) For 
the pp — >■ gd — >■ + X process, (b) For the pp gd t£f + X process, (c) For the 
— >■ —¥ tt~l + X process. The descriptions for the K-factor curves labeled (i), (ii), (iii) and 
(iv) are the same as in FigsJEJ 
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Figure 9: The NLO QCD correction to the LO cross sections (Aunlo) and the corresponding 
K-factors versus with A = — 1 TeV at the \fs = 7 TeV and \fs = 14 TeV LHC (a) 
For the pp — >■ gd — >■ + X process, (b) For the pp — >■ g<i — )■ + X process, (c) For the 
— >■ — >■ + X process. The descriptions for the K-factor curves labeled (i), (ii), (iii) and 
(iv) are the same as in FigsJHl 
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Figure 10: The NLO QCD correction to the LO cross sections (Aunlo) and the corresponding 
K-factors versus m~ tl with A = 1 TeV at the \fs = 7 TeV and \fs = 14 TeV LHC (a) For 
the pp — >■ gd — >■ + X process, (b) For the pp ^ gd t£f + X process, (c) For the 
— >■ —¥ tt~l + X process. The descriptions for the K-factor curves labeled (i), (ii), (iii) and 
(iv) are the same as in FigsJEJ 
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Figure 11: The NLO QCD correction to the LO cross sections (Aunlo) and the corresponding 
K-factors versus trig with A = TeV at the a/s = 7 TeV and -y/i = 14 TeV LHC (a) For 
the pp — >■ gd — >■ t£~ + X process, (b) For the pp — > gd — > i£f + X process, (c) For the 
pp — > gs — >• t£~ + X process. The descriptions for the K-factor curves labeled (i), (ii), (iii) and 
(iv) are the same as in FigsJHJ 
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Figure 12: The NLO QCD correction to the LO cross sections (Aunlo) and the corresponding 
K-factors versus irtg with A = — 1 TeV at the y/s = 7 TeV and y/s = 14 TeV LHC (a) 
For the pp — >• — >■ + X process, (b) For the pp — >■ — >• £i^~ + X process, (c) For the 
— > gs — > t£~ + X process. The descriptions for the K-factor curves labeled (i), (ii), (iii) and 
(iv) are the same as in FigsEJ 
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Figure 13: The NLO QCD correction to the LO cross sections (Aunlo) and the corresponding 
K-factors versus mg with A = 1 TeV at the y/s = 7 TeV and y/s = 14 TeV LHC (a) For 
the pp gd t£~ + X process, (b) For the pp — > gd — > t^/" + X process, (c) For the 
pp — > gs — > t£~ + X process. The descriptions for the K-factor curves labeled (i), (ii), (iii) and 
(iv) are the same as in FigsJHJ 
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Figure 14: The difference between the NLO QCD corrections with A = ±1 TeV and A = TeV/ 
versus at the = 7 TeV and = 14 TeV LHC. (a) For the pp ^ gd t£~ + X process, 
(b) For the pp — >■ — >■ + X process, (c) For the — >■ gs —¥ t£~ + X process. 
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Figure 15: The LO and NLO QCD corrected polarization asymmetries (A t ) of the (anti)top 
quark at the early and future LHC as the functions of mi with A = TeV. (a) For the 
process pp — )■ gd — )■ t£~ + X. (b) For the process pp — )• gd — >■ + X. (c) For the process 
pp —¥ gs —¥ t£~ + X. (d) For the process pp — > gs — >■ ti\ + X. 
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In further discussion, we consider the case where if = (i.e., i = 2). After the production 
of the scalar muon, there follows a subsequential decay of fi^ — > fi Xi with branch ratio 92.5% 
3l| . then the final state involves muon, the lightest neutralino and (anti)top-quark jet (/iXi^-(^))- 
As a demonstration we assume there exist = 189.9 GeV, A' 231 = 0.1 and the other A' = 0; 
then we present the LO, NLO QCD corrected transverse momentum distributions of the final 
(anti)top quark and the corresponding K-factors (K(p^) = da NLO {PT^) 1 1 d&LoiPj^)) f° r the 
processes of pp — >■ gd — >■ tjl~ — > + X aricl PP ~~ >" Qd — > — >■ iji + X\ + A at the LHC 

in Figs JTBT a.b.c.d). Figs flW a) and (c) are for the p^(Pt) distributions at the y/s = 14 TeV 
LHC for the processes pp — >• gd — >■ tf/,~Xi + an d pp ^ gd ^ + A, respectively, and 

FigsfTBTb) and (d) are at the -y/i = 7 TeV LHC for the processes pp — >■ — > tji^x°i + 
and pp gd ^ ifx + Xi + ^> respectively. These figures show that in FigsJT6l there exist peaks 
located at the position around p^(p^) ~ 80 GeV^ at the early and future LHC, separately. We 
can see that the LO differential cross sections are significantly enhanced by the QCD corrections 
with inclusive jet selection scheme, while the QCD correction by using an exclusive jet selection 
scheme keeps the convergence of the perturbative series in the plotted Pt(Pt) range. 

The transverse momentum distributions of the final muon and the lightest neutralino, and 
the corresponding K-factors for the pp gd — > tfi~Xi + X process at the \fs = 14 TeV LHC 
are depicted in Figs JT7T a.b). separately, while the corresponding distributions and K-factors at 
the a/s = 7 TeV LHC are depicted in FigslTTTcd). respectively. In FigsJT8Ta.b.c.d) we show 
the transverse momentum distributions and the corresponding K-factors of the final particles 
after the decay of jl + for the process pp — » gd — >■ + A. FigsJISTa) and (b) show the pr 

distributions of the final /x and Xi an d K-factors at the y/s = 14 TeV LHC, separately, and 
Figs JTBT c) and (d) demonstrate the pf^ and p^ distributions and K-factors at the yfs = 7 TeV 
LHC respectively. In the calculation for these results we use the narrow-width approximation 
(NWA) method to handle the resonant scalar muon effect. Here we assume A' 231 = 0.1 and the 
other A' = 0. The curves of K-factor labeled (I) and (II) correspond to adopting the inclusive 
and exclusive gluon/ ( ant i) quark jet event selection schemes, respectively. It is clear that with 
the exclusive jet event selection scheme we can keep the convergence of the perturbative series, 
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Figure 16: The LO, NLO QCD corrected differential cross sections of the transverse momentum 
of final (anti)top quark and the corresponding K-factors K{p®) = do jv lo {,Pt^ ) / dcr lo (p ) 
with A = TeV. (a) For the process pp ^ gd ^ t£~ + X at the y/s = 14 TeV/ LHC. (b) For 
the process pp — >■ grf — >■ t£^+X at the y^s = 7 TeV LHC. (c) For the process pp —> gd —> t£^~+X 
at the y/s = 14 Tel/ LHC. (d) For the process pp gd i£+ + X at the y/s = 7 TeV LHC. 
The curves of K-factor labeled (I) and (II) correspond to adopting the inclusive and exclusive 
scheme, respectively. 
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process 




^ = 50(3^7$? = 3 


^ = 35GeV^ = 2.7 


o"jvLo(fb) 


K-factor 


onlo (fb) 


K-factor 


pp — > gd 
-)• U~ + X 


501.0(1)/114.34(2) 


541.4(4)/134.9(2) 


1.081/1.180 


491.9(4)/123.6(1) 


0.982/1.081 


pp — > gd 

->• uj + x 


163.4(2)/28.39(4) 


184.7(4) /36.08(9) 


1.130/1.271 


169.7(3)/33.03(8) 


1.038/1.164 


pp — > gs 

«r + x 


97.78(1)/15.411(2) 


115.9(2)/19.85(9) 


1.185/1.288 


106.4(2)/18.07(8) 


1.089/1.172 



Table 3: The numerical results of the LO and NLO QCD corrected total cross sections by 
adopting the exclusive selection scheme for the pp — > gd — )■ U~ + X , pp — > — > t£ + + X and 
pp ^ gs tt~[ + X processes with A = TeV . The data on the left and right sides of the 
slash correspond to the results for y/s = 14 TeV and y/s = 7 TeV, respectively. There we take 
m- u = 189.9 GeV. 

and the NLO QCD corrections mostly enhance the LO differential cross sections of the final 
particles in the plotted pt range. 

In Table [3] we list some of the numerical results for the LO and NLO QCD corrected total 
cross sections by adopting the exclusive gluon/light-(anti) quark jet selection scheme for the 
pp —> gd —>• tt~l + X, pp — )• gd — )■ tt[ +X and pp — >• gs — >■ +X processes at the y/s = 14 TeV 
and 7 TeV LHC. There we take = 189.9 GeF, A- 31 = 0.1 and the other A' = for the 
first two processes, \' a2 = 0.1 and the other A' = for the last process. We consider the phase 
space with the restriction of either p^f* < Pj^j et = 35, 50 GeV or "if et > rfj£ = 2.7, 3 for 
gluon/light-(anti) quark jet transverse momentum. 

IV. Summary 

In this paper, we calculate the complete NLO QCD corrections to the single slepton associated 
with a (anti)top-quark production process at the early (y/s = 7 TeV) and future (y/s = 14 TeV) 
LHC. We investigate the dependence of the LO and the NLO QCD integrated cross sections on 
the factorization/renormalization energy scale, and study the influence of slepton, stop-quark 
and gluino masses on the NLO QCD corrected cross sections. We point out that the uncertainty 
of the LO cross section due to the introduced energy scale \i is apparently improved by including 
NLO QCD corrections, and the exclusive jet event selection scheme keeps the convergence of 
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Figure 17: The LO, NLO QCD corrected distributions of the transverse momenta of final 
muon and neutralino for the process pp — )■ gd — >■ t/i~Xi + and the corresponding K-factors 
(K(p T ) = with A = TeV at the LHC. (a) For the final muon at the yfs = 14 TeV 

LHC. (b) For the final x°i at the y/s = 14 TeV LHC. (c) For the final muon at the = 7 TeV 
LHC. (d) For the final x°i at the y/s = 7 TeV LHC. 
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Figure 18: The LO, NLO QCD corrected distributions of the transverse momenta and the 
corresponding K-factors (K(p T ) = da ^° /^tp-) of the final particles for the process pp — > 
gd ->■ tyu+x? + X with A = TeV at the LHC. (a) For final muon at the y/s = 14 GeV LHC. 
(b) For the final x? at the yfe = 14 GeV LHC. (c) For the final muon at the ^ = 7 GeV LHC. 
(d) For the final x°i at the ^ = 7 GeV LHC. 
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the perturbative calculations better than the inclusive scheme. We find also that the SUSY 
QCD correction generally increases if the t\ or g mass is getting large, but the non-zero tj—dk~U 
coupling with A = 1 TeV could distort that curve tendency in some cases. We present the 
LO and the QCD corrected distributions of the transverse momenta of final products involving 
(anti)top quark, muon and the lightest neutralino. Our results show that the NLO QCD 
corrections suppress the polarization asymmetries of final (anti)top quark, and mostly enhance 
the transverse momentum distributions of final particles. 
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